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A trajectory-preprocessing algorithm has been devised which matches an ten n a 
angular position , velocity, and acceleration to those of a target. This eliminates 
vibrations of the antenna structure caused by discontinuities in velocity and ac- 
celeration commands, and improves antenna-pointing performance by constraining 
antenna motion to a linear regime . The algorithm permits faster acquisition times 
and preserves antenna-tracking capability in situations where there would otherwise 
be an unacceptably sudden change in antenna velocity or acceleration. A simulation 
of DSS 13 shows that this preprocessor would reduce servo error to 1 mdeg during 
acquisition of a low-Earth-orbiting satellite. 


I. Introduction 

When a large antenna is moved in a sudden or jerky manner, the ensuing vibration of the structure 
can adversely affect pointing accuracy. A fast-moving antenna which stops in a precipitous manner upon 
acquiring a target may vibrate enough to lose lock. Even in situations for which the vibration amplitude 
does not immediately produce an unacceptable pointing error, the servo control may enter a nonlinear, 
and possibly unstable, region. This problem can arise when the servo controller attempts to apply an 
excessive velocity or acceleration in an attempt to track a sudden command change [1]. 

In the past, Deep Space Network pointing requirements have been for slow antenna motion. While 
DSN antennas are required t<^ track at up to 0.4 deg per sec, they are not required to meet pointing 
requirements for other than sidereal targets. 1 The prospect of very accurate acquisition of low-Earth- 
orbiting satellites (with angular velocities in excess of even 0.4 deg per sec) is a relatively new idea. The 
anticipation of such requirements has led to the present work. 

Trajectory preprocessing is one method used to preserve antenna-pointing integrity. The principal idea 
is to make antenna motion smoother during target acquisition, although the algorithm can be applied 
to all antenna motion commands. The present antenna control system is nonlinear. Its behavior is 
governed by the nonlinearity (acceleration and velocity limits), its inputs, the initial conditions, and 
the linear subsystem frequency response (the controller bandwidth). This has resulted in a system with 
complicated switching rules for changing the controller bandwidth as a function of inputs. A change in a 
velocity or acceleration limit affects the switching rules. The system can still work poorly for some sets 
of inputs, initial conditions, and parameters. 


W. Scherr, Deep Space Communications Complex Subsystem Functional Requirements , Antenna Mechanical Subsystem 
(1991 through 199 7 ), JPL D-1179, Rev. C (internal document), Jet Propulsion Laboratory, Pasadena, California, Septem- 
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Since commands which are in violation of these limits cause the system to behave nonlinearly, a simple 
method for allowing the system to operate in a linear regime is to force command angles to conform to 
the limits. Trajectory preprocessing performs this task. 

Even if an antenna were perfectly rigid and could be moved at will at accelerations of up to 1 deg/sec 2 , 
it could still overshoot a target during a high-speed acquisition. When the difference in speeds between 
average acquisition velocity and tracking velocity can easily exceed a deg/sec, some algorithm is needed 
to ensure that the acquisition time remains reasonable for a target which may not be visible for more 
than a minute or two in the first place. This provides a further motivation for considering trajectory 
preprocessing. 

The problem to be solved by this preprocessor is to find an optimal, or at least an adequate, path from 
an initial antenna position, velocity, and acceleration to some target trajectory. The first idea that was 
considered for finding this path was to use the calculus of variations. This idea was abandoned for three 
reasons. First, the method is overly complex. Second, the calculation time can be large. The method is 
iterative, and the calculation time is indeterminate, making it unsuitable for a real-time system. Finally, 
although the calculus of variations can give a least-time solution, it is not easy to include constraints that 
will prevent sudden changes in the slopes of the antenna velocity and acceleration profiles. 


II. The Three-Region Method 

The method described in this article involves three regions of antenna motion and was inspired by the 
following scenario. A target is far away and one wishes to move the antenna towards the target trajectory 
as quickly as possible. So one begins by accelerating the antenna to its maximum speed; that is region 1. 
Then, in region 2 one moves the antenna at maximum speed until one is near the target. Finally, in 
region 3, the antenna is decelerated until it matches the apparent target angular position and angular 
velocity. 

As applied to a low-Earth-orbiting satellite, the initial scenario envisioned the following preprocessor 
steps: 

(1) Input a set of target positions and velocities (9 and v) for both local azimuth and elevation as 
a function of time using known orbital parameters (updated by optical detection). 

(2) Choose an intercept (acquisition) 9 e/ = 9 0cl . From this obtain the time interval, T, from the 
start of region 1 to the intercept time, as well as 6 az and both v az and v e j. Input the initial 
antenna pointing (do and t>o) nnd the maximum angular velocity for the antenna (iVnax)- Input 
or calculate maximum angular accelerations for the antenna (a, max > 0). 

(3) If necessary, try a different maximum acceleration (remaining less than or equal to the specified 
limit) or a later acquisition time. Set a flag if acquisition is not possible before 6 e \ > 0 max . 

The same algorithm can be used for reacquisition during tracking (without the flag). One does not 
need to start by choosing the acquisition elevation. A total acquisition time can be input instead. This 
time can be chosen by noting the distance in position and the change in velocity that must occur. Using 
one’s experience from previous acquisitions, it will often be possible to choose an acquisition time within 
a second of the minimum. Should one’s choice be too small, the algorithm can be rerun with steadily 
increasing total times input. In practice, the maximum acceleration may simply be set to the specified 
acceleration limit. 

For acquisition, in both azimuth and elevation, it is best to chase the target rather than approach it 
head-on; this lessens the required change in velocity. That means that 0q should be a little less than 0. 
Thus, while waiting for the target to appear, the antenna will generally be pointing at a place in the 
trajectory close to that anticipated by an optical acquisition aid. The antenna will typically have an 
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initial velocity of approximately zero in elevation (or it might badly overshoot the target). When the 
target shows up “late,” the Earth’s rotation will cause it to move in azimuth as a function of lateness, so 
the antenna may initially be moving slowly in azimuth. The formula for the minimum time to acquire is 
then T — vj/A y where A is the average acceleration. 

When A = a max , we get the fastest acquisition. An a max in elevation of 0.5 deg/sec gives T = 2 vj, 
where v is in deg/sec and T is in sec. 

In our acquisition schemes, the average A will often be a max /2 or less, giving a minimum acquisition 
time of 4vj sec. When vj is about 0.5 deg/sec, this results in a minimum acquisition time on the order 
of 2 sec. 

This minimum acquisition time applies only when 


v f 

A ° el = 2 A 


If, for example, we initially aim at the acquisition point instead of a point Vj/2A below it, the minimum 
time increases from T to Tj\f 2 — 1. 

Figure 1 shows the matching in elevation of antenna pointing with that of a target. Figure 2 shows 
the antenna angular velocity for this example. In these figures, the antenna is aimed a little above 
the acquisition point. Figure 3 shows how the algorithm matches angular position when the antenna is 
initially aimed a little below this point. 

The acquisition time will often be determined by the change in v e {. The final azimuthal velocity is 
usually much smaller, and there may even be an initial azimuthal velocity. However, if the antenna is 
simply waiting several degrees from the acquisition, 0 az , the azimuthal acquisition time can easily become 
the total acquisition time. 


III. Acquisition Schemes 


Four acquisition schemes are considered in this article: 

Scheme 1: Match initial and final angular position and velocities. Use the maximum 
acceleration throughout each acceleration region. 

Scheme 2: The same as scheme 1, but use a sinusoidal acceleration pattern to avoid 
continuities in acceleration. 

Scheme 3: The same as scheme 2, but match the final acceleration as well. 

Scheme 4: The same as scheme 2, but match both initial and final accelerations. 


allowable 
large dis- 


Each of these ideas has some merit. The calculations needed to implement the trajectory-preprocessing 
algorithm are similar for any of the first three schemes. 

A. The First Acquisition Scheme 

Scheme 1 has the advantage of speed. By using the maximum acceleration, the target can be reached 
more quickly. Even when the increased initial overshoot and ensuing oscillations are taken into account, 
the total acquisition time may be minimized in some situations. On the other hand, when the angu- 
lar distance to the target is large and the maximum allowable pointing error is small, this scheme is 
inappropriate. 
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Fig. 1. Matching of antenna and target positions for 
constant accelerations (elevation only). 



TIME, sec 

Fig. 2. Matching of antenna and target velocities for 
constant accelerations (elevation only). 



0 2 4 6 8 10 12 14 16 

TIME, sec 


Fig. 3. Matching of antenna and target positions for 
elevation only (chasing target). 


This scheme uses the maximum acceleration in each acceleration region. Although this introduces 
discontinuities in acceleration, it permits a straightforward calculation of the acquisition parameters. 
Then a = ±a m = ±a mar in regions 1 and 3. By integrating twice, we get the velocities and positions as 
a function of time in all three regions. The following calculations must be performed both for elevation 
and for azimuth: 


For region 1, 


v = t>o + a 1 1 


( 1 ) 


For region 2, 


9 — 6o + vot-\- 


a\t 2 

~ 


v — V2 = vo + a\t\ =vj— a 3 ts 


( 2 ) 

( 3 ) 


9 = Oq + V 0 ti + X)2 (t — fl) + 


a\t 2 i 

2 


( 4 ) 
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where T = *i + t 2 + t 3 . 

For region 3 , 

v = vj - a 3 (T - t) ( 5 ) 

(T - t ) 2 

e = e f -v f (T-t) + a 3 K -^-L (6) 

1 . Calculation of Acquisition Parameters. In the above equations, we do not yet know the signs 
of a\ or a 3 . Nor do we know the durations, or f 3 , or the constant velocity, v 2 . These are calculated as 
follows. Input the initial and final antenna positions (azimuth and elevation), Oq and Of, 

AO — Of — 0 o 

as well as the initial and final antenna velocities (azimuth and elevation), vq and v j , 

Av ~ Vj — Vq 


and the total time, 


T — - 1 1 in 1 3 


as well as the maximum acceleration, a m , and velocity, v m (in both azimuth and elevation), which must 
be less than or equal to the requirements, a max and v max (in both azimuth and elevation). 

From the above inputs, output and u 2 and also determine a\ and 03. To obtain these outputs, 

the following formulas are used. First, a normalized angular position, x, and a normalized angular velocity, 
y, are calculated: 


AO uo 

a m T 2 a m T 


( 7 ) 


y = Av 
a m T 

Next, and £j, intermediate variables that are used to find the signs of the accelerations in regions 1 
and 3 , are determined: 


V 2 V 2 

£0 =£/=—! when y < 0 and y + — < x < — — 


y 2 y 2 

£0 = £/ = 1 when y < 0 and — < x < —y — — 


y 2 y 2 

£o=l and £j — — 1 when y > 0 and x > y — — or y < 0 an< ^ x > 


y 2 y 2 

£0 = — 1 and £/ = 1 when y > 0 and x < — or y < 0 and x < y + 


( 9 ) 

( 10 ) 

(11) 

( 12 ) 
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These give us a i and a 3 : 


— £0^m 

(13) 


(14) 


Now the velocity in region 2 can be found: 


fore 0 = £/, V 2 


y 2 g/ - 2z 
2 (ye f - l) 


for eo # £/ , 


y2 


ye/ - 1 + yVepe/ - 2 ye/ + 2x(c; - e 0 ) + 1 
€j —€ 0 


U 2 = + V 0 


Finally, the time intervals can be deduced: 


v 2 - vo 

f l = 

a i 

<3 

^ 2 = T ~ f 1 — <3 


(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 


Once the above parameters have been obtained, position and velocity commands are determined from 
Eqs. (1) through (6). 

The main problem with this acquisition scheme is the discontinuity in antenna acceleration at the 
borders of regions 1 and 3. The problem is most serious at the border of region 3, where transient 
phenomena may significantly affect antenna-pointing accuracy immediately after acquisition. To see this, 
consider the example of Figs. 1 and 2. In this, Example (1), 

9 0 = 25.104 deg 
0 f = 24.253 deg 
i? 0 = —0.001 deg/sec 
Vf — 0.479 deg/sec 
T — 6.6 sec 
\a\ = 0.25 deg/sec 2 

which has the solution 
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V 2 = —0.46 deg/sec 
<1 = 1.85 sec 
^ 2=1 sec 
t 3 = 3.75 sec 

Figure 4 shows anticipated servo error during and immediately following acquisition for this example 
using a simulation in which a PI controller is applied to a model of a 34-m Deep Space Network antenna, 
namely DSS 13 in Goldstone, California. This servo controller, which uses both proportional (P) and 
integral (I) feedback terms, is in use at DSS 13. The controller and the model of the antenna are described 
in [ 2 ]. 



-6 1 1 — I — i uJ , , , * I 

0 2 4 6 8 10 12 14 16 

TIME, sec 

Fig. 4. Anticipated servo error for Example (1) with 
constant acceleration in regions 1 and 3 (elevation 
only). 


2 . Derivation of Results. This section gives a derivation of the results of Eqs. (7) through (20). 
Begin by conserving the total angular distance: 


V\ti + 4" ^3^3 — A# 


( 21 ) 


The values t\ and t 3 are obtained from the change in velocity divided by the acceleration. This immedi- 
ately gives Eqs. (18) and (19). 

The values v\ and U 3 are the average angular velocities in regions 1 and 3. Since the acceleration is 
constant in each of these regions, 


vi 


V2 + Vp 

2 


and v 3 = 


V2 + Vf 

2 


When Eqs. (18), (19), and (22) are substituted in Eq. (21), the result is 


- ^0 

2 a, 


+ + 



2a 3 


= AS 


( 22 ) 


(23) 
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By substituting T - t x - t 3 for t 2 in Eq. (23), we get 


_ (v 3 - i’o) 2 _ + 2nT + (V ~ = 2Ae 

a\ a 3 


(24) 


Further substitution of the definitions of x,y, 1 / 2 5 ^ o > and £/ from Eqs. (7), (8), (17), (13), and (14) gives 

yl{ej - £ 0 ) + 2 y 2 (l - ytj) + 2 / 2 ^/ - 2x = 0 (25) 

When £0 = £/, which is often the case, Eq. (15) follows directly from Eq. (25). When £ 0 i 1 £/. Ec l* ( 16 ) 
immediately results from applying the quadratic formula to Eq. (25). However, the sign in front of the 
square root in Eq. (16) is yet to be determined. This is done as follows. Note that 


h = T- 


a 3 a\ ai a 3 a 3 


vo_ 

a 3 


which gives 


~ = 1 -ysf+ yiiej - £ 0 ) 

The square root term in Eq. (16) is also equal to 1 — y£j -j- 2 / 2 (^/ ~ £ 0 )* Since ti/T must be positive, the 
plus sign must be chosen in front of the square root. 

The next step is to find out where the solutions for y 2 described in Eqs. (15) and (16) are valid. The 
values of x and y for which a solution can be found will be referred to as an “area of validity,” which 
can be plotted on a graph of y versus x . This area is in ‘phase space, and as it increases, a larger 
number of combinations of target positions and velocities can be matched. Starting from Eq. (16), the 
minimum and maximum values of £ as a function of y 2 are obtained by setting dx/dy 2 = 0. This gives 
y 2 = (j, £/ _ i)/( £/ _ £o ) and t 2 = 0. This means that the discriminant of the square root in Eq. (16), t\, 
is 0 at the borders of the area of validity of solutions to Eq. (25). Setting the discriminant to zero gives 

_ 1 - 2 yej ± y^£oef (26 ) 

2(e 0 -£/) 


At the external boundaries of the area of validity, e 0 f £j . Thus, the area of valid solutions is bounded 
by the curves 


y 2 + 2y - 1 


(-l<2/< 1) 


and 


I= l + 2p! (-,<„<!) 

This area of validity is shown in Fig. 5. Notice that solutions for which x and y have the same sign are 
favored (have more phase space) than those which do not. 


146 



Fig. 5. Area of valid solutions ( a m = s max /2 
and a Q = a^= 0). 


In these calculations, we have ignored the situation in which v max corresponds to |t/ max | < 1, which 
would cut off the region of valid solutions aty< 1 or y > —1. This can be dealt with as follows: Note 
that for the antenna velocity to exceed a velocity limit, either v 0) uy, or v-i must exceed that limit. If 
t>o exceeds the limit, we already have a problem, but it is one that has little to do with the trajectory 
preprocessor. If vj exceeds the limit, we cannot acquire at that point, but we can try other values of T 
for which vj may be smaller. If v*i is the only culprit, we increase T to a point at which t>2 is acceptable. 

Figure 5 also shows internal borders in the area of validity. These borders define the regions where 
a\ and are both negative, both positive, of opposite sign with a\ positive, and of opposite sign with 
a i negative, respectively. At these borders, t\ — 0 (which gives 2/2 = 0) or = 0 (which gives y 2 = y). 

Substituting in Eq. (25), 


y 7 = 0=>X= 2 / 

(27) 

. y 2 e 0 

j/ 2 = j/=>I = y 

(28) 


By defining the borders of each of the four internal regions, Eqs. (27) and (28) permit us to write down 
Eqs. (9) through (12), completing our derivation of the formulas used in the first acquisition scheme. 

B. The Second Acquisition Scheme 

To avoid the discontinuities in acceleration, one can choose a sinusoidal acceleration pattern. This 
trades acquisition rate for pointing accuracy; the price that is paid for improved accuracy is a factor of 
two in average antenna angular acceleration. The preprocessor matches position and velocity as before. 
The antenna angular acceleration is still not matched to that of the target, but it is zero rather than the 
maximum allowable acceleration at acquisition. The idea is to let a — ±a m (l — cos2ttu;) for some u? and 
for a m < a max / 2. In particular, for region 1, 


Integrating with respect to t, 


a. — o,\ 


2ttA 

— cos ] 

1 1 ) 
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(29) 


Integrating again, 


v — v Q -\- a 1 



2 x 


sin 


2xA 

w 


6 = O 0 + v 0 t + a\ 



t\ tj 2wt\ 

j? + 4 ^ cos ir) 


(30) 


For region 2, 


a — 0 


v = V 2 = vq + ait i = vj — a 3 T 3 


(31) 


0 = #0 + Voti + V 2 (t — 1 1 ) + 


ait\ 

2 


(32) 


For region 3, 


a = a 3 


^1 — cos 


2tt(T-0\ 

^3 ; 


Integrating with respect to (T — i), 


v — vj — a 3 



U . 27t(T - <)' 

— — sin 

27T ^3 , 


(33) 


9 = 9j — vj(T — f) + <33 


tj t\ 

\ 2 47t 2 47t 2 


27r(T-<)\ 
<3 ) 


(34) 


Note that da/dt is not merely finite everywhere; it goes to 0 at the borders of each region. 

At first, these equations seem more complicated than those for constant accelerations, but the solutions 
for their parameters are identical. Equations (7) through (20) still hold. The acceleration is not constant, 
but Eq. (22) is still valid. The maximum value of a m is reduced by a factor of 2, but Eqs. (13) and (14) 
are unchanged. The same four solutions exist with the same boundaries in x and y, as represented by 
Fig. 5. The only difference is that Eqs. (1) through (6) have been replaced by Eqs. (29) through (34). 

Example (1) still has the same solution for t x , t 2 , t 3) and V 2 . However, the anticipated servo error has 
been reduced because the acceleration changes more smoothly. 

Figure 6 shows the matching of antenna and target angular positions in Example (1) for acquisition 
scheme 2. Figure 7 shows the corresponding matching of velocities. Note that the velocity slope is zero 
at acquisition for the antenna but not for the target; there is still a slight discontinuity in acceleration 
at acquisition. Figure 8 gives the anticipated servo error for Example (1) using the second acquisition 
scheme. 
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Fig. 6. Matching the target position for raised cosine 
acceleration (elevation only). 


Fig. 7. Matching the target velocity for raised 
cosine acceleration (elevation only). 



Fig. 8. Anticipated servo error (for Example (1) with 
raised cosine acceleration in regions 1 and 2 (elevation 
only). 


Since T has already been chosen prior to the calculation of t\ and £3, it may seem unnecessary to set 
a m to a max /2 when a smaller a m would suffice. Also, potentially huge values of da/dt (jerk) can arise 
in this trajectory preprocessing scheme. These high values of jerk are produced when one is near the 
internal borders within the region of valid solutions in Fig. 5, since <1 or ts approaches 0. Since | a | 
goes from 0 to 2a m and back in each region, the average \da/dt\ is Aa m /ts in region 3. The maximum 
\da/dt\ is 8 a m /t 3 in region 3, and as ts approaches 0, this number becomes larger until the increase is 
filtered out by the servo mechanism. Actually, this gives us a very small servo error because for the PI 
controller, the maximum error caused by a sudden pulse does not exceed the total displacement caused 
by the pulse. Thus, the PI controller does not react much to a ts of 0.1 sec because 1/2 a av t 2 for an 
average a of 0.125 deg/sec 2 is only 0.6 mdeg, which is an upper bound on the servo transient error. This 
is illustrated by Example (2), for which 

0 O = 23.7618 deg 
6 j = 24.253 deg 
vj — 0.479 deg/sec 
v 0 = —0.001 deg/sec 
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H = 

0.25 deg/sec 2 

T = 

2 sec 

<i = 

1.9 sec 

t‘l ~ 

0.0 sec 

^3 = 

0.1 sec 

V2 = 

0.4665 deg/sec 

Figure 9 shows 

the anticipated 



Fig. 9. Anticipated servo error during acquisition for 
Example (2), (f 3 = 0.1 sec, f 2 = °) (elevation only). 

If it is thought that the servo controller would react adversely to such a sharp acceleration pulse, one 
can modify d m so that neither t\ nor £3 is small, or at least choose a reasonable a m to begin with and 
modify either a m or t if necessary. For example, one can choose an a m that gives t 2 = 0 and if ti or t 3 is 
small anyway, one simply increases T. From Eq. (26), t 2 = 0 gives 


Picking £ 0 £/ = — 1, 


Let 


Then 


y 2 e 0 £j - 2ye; + 1 = 2 (e 0 - £j)x 


y 2 + 2ye j f 1 “ j) ~ 1 = 0 

2x A vT - 2A 9 + 2 v 0 T 
~ 1 “ V _ A vT 

-% = v = -ej n±VwTi 


(35) 


The sign in front of the square root in Eq. (35) and the value of Sj are readily determined by observing 
that y has the same sign as Ais and that |r/| < 1. This determines y and gives 

Au 

a rn — ~77V 

yT 
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The a — 0 region can also be omitted from the algorithm entirely. The resulting scheme would have 
no intermediate region of maximum velocity for large changes in antenna position and would increase 
acquisition time as a result. If by some chance, 

A y _ vj - vl 
~¥ ~ 2A 0 


(and even with some latitude in choosing t f it is unlikely that this will occur), then a one-region solution 
is possible where 


a — 


2L 

Av 


27 rt 
~T~ 


) 


If not, then SqSj = — 1, and y is determined from Eq. (35), which also gives a m} while 


By leaving out the a = 0 region, overall accelerations are reduced. However, there is an increased flexibility 
in maintaining a three-region algorithm, and picking a m prior to solving for j /2 may be impractical in 
some situations. 

C. The Third Acquisition Scheme 

In our examples, the target acceleration and servo controller are such that the sudden discontinuity 
in acceleration at acquisition does not significantly increase the pointing error. However, there is a 
straightforward way to avoid this discontinuity by modifying the previous acquisition scheme so that it 
matches the target acceleration at acquisition instead of acquiring with zero acceleration. One simply 
switches to a frame of reference which has an acceleration equal to that of the target. Figure 10 shows 
the anticipated servo error for Example (1) in this case. 



0 2 4 6 8 10 12 14 16 


TIME, sec 

Fig. 10. Servo error for acquisition scheme 3 and 
constant velocity track. Acquisition is at t = 6.6 sec. 

Let ay be the acceleration of the target at t = T. Then match a position of 9 } - 1/2 aT 2 and a velocity 
of Vf — aT. 

This means that Ad now equals 0/ - 0 O - ((1/2 )aT 2 ) and Av - vj - v 0 - aT. ’Once again, Eqs. (7) 
through (20) remain valid. However, care must be taken to ensure that the actual antenna accelerations 
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and velocities never exceed the maxima. The acceleration can be kept within bounds by choosing a m = 
a max - |a/ 1/2. But the velocity is trickier, as the actual velocity in region 2 is no longer constant. For a 
situation in which the antenna must traverse a large angular distance to acquire a target, it might at first 
appear that the maximum velocity will be reached at the end of region 2. This is not the case. At the 
end of region 2, the antenna has an acceleration of aj , so the velocity is still increasing. The maximum 
velocity is reached when the acceleration first reaches zero in region 3. 

The difference between this scheme and the first scheme is that Eqs. (1) through (6) are now replaced 
with the following equations: 

For region 1, 


a — a,\ 


(. 2ttA 

^_ cos __J +a/ 


. 1 1 . 2 7rt 

v = Vo + ai(t — — sin h a ft 

27T fi 


(36) 


6 = 0 O + vot + ai 



/ 2 t 2 
l \ , 1 i 

47T 2 47 r 2 


2 7 Tt\ 
COS I 

h ) 




2 


(37) 


For region 2, 


a — a j 


v = i’2 T a jt 


(38) 


0 = 0q 4- v Q ti + — ^i) T 


2 


+ 


ayt 2 

~T” 


(39) 


For region 3, 


0 — 0j - Vj(T -t) + a 3 


cos 

2tt(T — f)^ 

*3 ) 

f 

T - 

. h 

- 1 — — sin 

\ 

27T 

(T- 

-0 2 4 


+ U/ 

2t r(T-ty 

^3 




<3 

4^2 + 4^ COS <3 


2t v(T-t)\ t a,* 2 
1 


+ 


(40) 

(41) 


D. The Fourth Acquisition Scheme 

Our final scheme applies to situations where the antenna has an initial acceleration ao and a final 
acceleration a;. In this scheme, both the initial and final accelerations are matched by the preprocessor. 
We will assume that an initial value of a m is chosen either by setting it to a max /2, by looking at Eq. (35) 
(possibly modified by adding or subtracting a 0 or a;), or by some as yet undetermined method. Inputs 
are then 0 Ql v ai v f> A 0 , Aj } A m , T, v maX} and a max (for both azimuth and elevation). Equations 
(1) through (6) are now replaced by the following: 
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For region 1 , 


a — a\ [ 1 — cos 


2ttV 


«rM( ,+ ~?) 

Note that da/dt = 0 both at t = 0 and £ = ti. By integrating, 

v = vo+(a 1 + ^y + ±- (a 0 sin y - a, sin 

= 0o + v 0 1 + («i + y) y + [2a 0 (l - cos y) - <»i (l - cos ^ 


0 = 


For region 2, 


a = 0 


For region 3, 


u — ^2 - Vo + + y) t\ = v f - ^a 3 -j- y) *3 

0 = 00 + t>o*i + «,(* - <0 + («1 + y) f + 


/\ 2w(T — t)\ a/ 

a = a 3 ( 1 - cos — ] + ~4- 

h 




(42) 

(43) 


(44) 


(45) 


(„ , a S\ (T A h ( . *(T-t) . 2Tv(T-t)\ 

v ~ V J ~ ( a 3 + yj (T - t) - — sm ^ - a 3 sm — ^ 1 


(46) 


0 = 0 , 


- V ,(T - t) - (« 3 + ± [ 2o/ (, _ cos ££-Li>) - „3 (l - cos 


(47) 


When sign (a 0 ) = sign (ai), 


|nj j — ^ (2 a m |a 0 |) 

Otherwise |ai| = a m . When sign (ay) = sign (a 3 ), 

l a 3| = ^ (2a m - |ay|) 


Otherwise |a 3 | = a m . 
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I III) Ill ' ll III II I 


The solution from the previous acquisition schemes no longer applies. Although y 2 has a similar form, 
e 0 and tj are changed. The parameters are derived as follows. Once again, v\t\ + v 2 t2 + ^3 = A0. Let 

1 

Po ~ a l+ (a 0 /2) 

1 

0} ~ a 3 + (ay/2) 

Then 


t i = /3q{v2 ~ v 0 ) 
<3 = /3y(vy - d 2 ) 


From Eq. (45) we can calculate the change in position from t = 0 to < - h. Dividing by 1 1 , we get »i: 


v\ 


Vq + ^2 Mo 
^ r •? 


Equation (22) is no longer valid, so the calculation becomes more complex: 


V3 


vj + v 2 haj 


Plugging in, 


{v f - u 2 ) 


2 


/?/ + 



- (d 0 - D 2 ) 


2 


2#fao] 

it 2 J 


= 2A0 - 2v 2 T 


Let 


Cj — ci m 



a — (ft _ 
and Cq — flm l kO ^2 I 


Then 


(Dy - d 2 ) 2 e y - (d 0 - d 2 ) 2 c 0 = 2a m (Ad - v 2 T) 


Changing variables, 

y 2 (ey - Co) - 2y 2 (ycy - 1) + y 2 cy - 2z = 0 


When cy = c 0 , we get Eq. (15). Otherwise, 
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y 2(<7 - Co) = yej - 1 ± yjy^e 0 e f - 2 ye/ + 2a;(e / - e 0 ) + 1 


(48) 


which is the same as Eq. (16) except that the sign in front of the square root can be negative. 

Since eo can have at most 2 values, €j can have at most 2 values, and the sign in front of the square 
root can have at most 2 values, at worst one needs to solve Eq. (48) eight times to see if a valid value of 
2/2 can be obtained. In practice one does not need to try all eight possibilities. If y > 0, and cj cannot 
both be negative. If y < 0,Co and cj cannot both be positive. So only six cases remain. 

It is tempting to simply write down the equations for the borders of each region and pick which solution 
is valid. That would determine to and cj. However, this method is not practical in general, as will be 
illustrated by trying it for a 0 = — a m and aj = a m . When a\ and 03 are both negative, 

2 8 

eo = -1 + —=■ and e, = -2 + — 

7 t z 7 r J 


The minimum x is found by setting to y. From Eq. (48), 


x = y - 


o 

2 


Here, — 1 < y < 0. When —1 +ct//2a m = — 1/2 < y < 0, the maximum x is determined by setting y 2 — 0. 
Now Eq. (48) gives 



For -1 ^ y < — 1/2 — — l-|'Q//2a m , the maximum x is found by setting y 2 — l-j-2 y. Here Eq. (48) gives 

* = y 2 + f (<7 + 6 ) + ^( 3« 0 + 4 ) 

When a 1 and a 3 are both positive, 

8 2 

Co — 2 H — - and €t = 1 4 — - 

7 T 2 7 T 2 

Now it takes three curves to describe the minimum x. For 

y 2 

then 7/2 = 0 => x = —cj 

y 2 

then y 2 = y => x = y - — e 0 

then j/2 = 1 - y => x = y 2 {2tj - 3)y + ^(4 - 3f;) 


0 < y < 


Co + Cy 


2 do 

— t — < y < -h 1 T - — , 
Co + c/ 2a m 


1 + -p- <y< i, 

4 Clm 
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The maximum x also requires three curves. When 


0 < y < — , then 
co 


y 2 - y ^ X - y - 


y 2 e o 


1 1 . 

— < y < — , then 
£o € J 


ye/ ~ 1 . y 2 c of/ - 2 yej + 1 

y2 = T7^^ a:= 2(fo — f/) 


— < j/ < 1, then 
<7 


2/2 = 0 =S> X 


_ y 2< / 


When «i is positive but a 3 is negative, the minimum and maximum x require three curves each. 


8 8 

e 0 = 2 + — and e; = -2 + — 

7T 2 


The minimum x is as follows. For 


- r < v < 1 + Q ° . then 

2(e 0 - 1) 2 a m ’ 


2y+l l-2ye/-4y 2 

* =^-=- I = 8 


0 < y < 


2(eo - 1) 


, then y 2 = y => x = y - 


«/ 2 fo 


-1 + TpC- < y < 0, then 


y 2 = 0 =» X = 


y 2 c/ 


The maximum x is as follows. For 


1 ao . 

— < y < 1 + - — , then 
Co 


y 2 = y => x-y- 


V 2 e o 


0 < y < — , then 
Co 


l -yet l -2ycj + ej€ 0 y 2 

2/2 = 1 =>* = 3 


-1 + < y < 0, then 

2 dm 


V2 = 


2y + 1 1 - 2 ye/ - 4y 2 


8 


When a\ is negative but a 3 is positive, we obtain our final six curves. 


2 2 

fo = — 1 H — o and € J ~ 1 + — 2 

TT " ')T“ 


The minimum x is as follows. For 
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1 yrtf 

— < y < 1, then 2/2 = 0,=J> x - — — 
( J 2 


0 < y < — , then 
€ f 


2/2 


ycj - 1 _ -J/ 2 e/fo + 2 ytj - 1 

- => x — 

2 4 


— 1 < y < 0, then 


2/2 


y - l . y 2 + lyt} - i 

=>• I = 

2 4 


The maximum x is as follows. For 


< y < 1, then 

2e/ — 1 


2/2 


2 / — 1 y 2 + 2yc f - 1 

=> X == T 

2 4 


0 < y < 


2c j - r 


then 


2/2 = 0 =» ^ = 


y 2 tj 


— 1 < y < 0, then 


y 2 ^o 

2/2 -y^x = y — 


For other values of ao and a/, different curves determine the boundaries of the regions. The difficulties 
in finding the internal borders of the area of validity are sufficiently great that were this algorithm to be 
implemented , one would simply try all six candidate solutions and pick the first one that worked. 

Figure 11 shows the area of valid solutions for a 0 = — a m and aj — a m superimposed on the one for 
a 0 = af =0 shown in Fig. 5. Note the loss in phase space for valid solutions in the aj = — ao — a m case. 



Fig. 11. Shrinkage of region of valid solutions for 
a f = = a m- 


IV. Discussion and Conclusions 

Acquisition scheme 4 is overly complex even if one does not calculate the internal borders of the area 
of validity. Such an algorithm would be very difficult to implement and maintain. Scheme 3 should 
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suffice for high-speed acquisitions and is recommended as an option on Deep Space Network antennas. 
The other schemes produce acceleration discontinuities which may cause undesirable excitations of the 
antenna structure. 

The algorithm can be run with the entire set of commands output at once. However, in practice, it 
is not necessary to calculate all commands before implementing the first one. It is sufficient to find an 
acceptable acquisition time, T. The position and velocity commands can be calculated in real time as 
the antenna moves. If it is inconvenient to calculate the commands in real time, there may be enough 
processor time to calculate them in advance and store them. 

Since the trajectory preprocessor cannot be used to supply commands until a satisfactory acquisition 
time is calculated, the question arises of how to find the acquisition time quickly. The problem is not 
trivial, as the target may have a trajectory that is very difficult to match. There are a number of 
possible strategies for picking a candidate acquisition time. A trade-off is involved. If the processor is 
so slow that it may take several tenths of a second to discover if a candidate acquisition time will work, 
one must be conservative in one’s choice of candidate acquisition times. Only a few candidates can be 
tried. The candidate which is finally selected may not be optimal, but the time lost in finding a better 
solution may more than make up for the time saved by the improved answer. On the other hand, if 
thousands of candidate solutions can be tried in a second, one should expect to find an acquisition time 
that is within a fraction of a second of optimal. Mere processor speed does not guarantee success, as the 
processor may need to be shared with other tasks. Nor should it be forgotten that each target position 
may need to be translated in a relatively time-consuming manner to an equivalent antenna command. 
There may be coordinate conversions to apply, refraction must be taken into account, and subreflector 
squint must be corrected for, along with a host of other tabled or modeled systematic pointing errors. 
A trajectory preprocessor algorithm should not be developed for an antenna system unless one has a 
reasonable knowledge of the required and available processor time, both before and during the time 
period in which preprocessor commands are to be output. 

When the trajectory preprocessor is used to match the trajectory of a sidereal object, the acquisition 
time is easy to estimate. There is usually no hurry to acquire the object, but even if there were, the 
preprocessor could handle the situation rather easily. A problem arises when the algorithm is used to 
match the trajectory of a fast-moving object that may be moving at 10 to 50 percent of the maximum 
antenna angular rate at acquisition. 

It is understandable that one might wish to construct a simple, general algorithm to generate candidate 
acquisition times, One can guess a time of 1 sec, and should that be insufficient, continue with guesses 
of 2, 4, 8, 16, 32, 64, and 128 sec until a solution is found. If powers of two seem inappropriate, one 
can try powers of 3, 1.5, 1.2, or whatever. One can be satisfied with the first acceptable solution, or 
one can backtrack, looking for an even better one. Another idea is to start with the maximum A 0jv max 
and 2Au/a 7Tiair . Any of these ideas may be acceptable, but it seems far better to produce a carefully 
constructed table of candidate acquisition times for each given mode, especially for low-Earth-orbiting 
satellites. A simple default mode can be included as well. 

This trajectory preprocessor algorithm depends greatly on the ability to predict exactly where the 
antenna will be when preprocessor commands are to begin. Any confusion resulting from misapplication 
of pointing corrections or differences between hoped for and actual position will result in a discontinuity 
in command position, which is precisely the problem that trajectory preprocessing is supposed to avoid. 
Choosing anticipated commands for the initial antenna position and velocity does no good if the antenna 
is pointed elsewhere. Using the actual antenna position does little good if the antenna is moving quickly 
and the preprocessor commands are due to start only a second or two later. Other initialization errors 
are possible that could render the preprocessor ineffective. For example, if incorrect or inappropriate 
velocity or acceleration limits are used, so that the antenna cannot respond properly to the preprocessor 
commands, antenna control will be back into the nonlinear region that the preprocessor was designed to 
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rescue it from. As long as care is taken to avoid mistakes of this sort, the preprocessor should serve a 
useful function. 

There is no reason to demand that an antenna controller always be in a mode for which trajectory 
preprocessing is in use. The preprocessor can be an option used especially in cases for which very accurate 
tracking is required or for which control problems are anticipated. For this reason, an acquisition scheme 
that matches target acceleration may be favored over one that does not; the preprocessor may be used 
primarily when one wishes to avoid what sometimes seem like small acceleration discontinuities. However, 
use of a preprocessor as an option does not mean that it should be considered as an ad hoc feature rather 
than an integral part of a control system design. The issues of how and where to fit preprocessing into 
a system should be addressed even if it is not yet decided whether or not such an algorithm will be 
implemented. This would avoid problems that may arise when one attempts to add it after the rest of 
the system is complete. 

As pointing requirements become more strict and tracking speeds increase, trajectory preprocessing 
will become a more and more valuable option to improve antenna control. The algorithm described in 
this article could be put to good use in the Deep Space Network. 
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